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We derive a constraint on patchy screening of the cosmic microwave background from inhomoge- 
neous reionization, using off-diagonal TB and TT correlations in WMAP-7 temperature/polarization 
data. We interpret this as a constraint on the rms optical-depth fluctuation Ar as a function of 
a coherence multipole Lc- We relate these parameters to a comoving coherence scale, of bubble 
size Rc, in a phenomenological model where reionization is instantaneous but occurs on a crinkly 
surface, and also to the bubble size in a model of "Swiss cheese" reionization where bubbles of fixed 
size are spread over some range of redshifts. The current WMAP data are still too weak, by several 
orders of magnitude, to constrain reasonable models, but forthcoming Planck and future EPIC data 
should begin to approach interesting regimes of parameter space. We also present constraints on 
the parameter space imposed by the recent results from the EDGES experiment. 



Introduction — Simple estimates have long shown that 
radiation from the first star-forming galaxies in cold- 
dark-matter models should reionize the intergalactic 
medium (IGM) at redshifts z ^ 10 T, but the details of 
this epoch of reionization (EoR) are still unclear. Quasar 
observations suggest that the tail end of reionization oc- 
curred at a redshift ^: ^ 6 [2], but the implications of 
these measurements are difficult to interpret precisely ^3J . 
A constraint r = 0.074 ±0.034 to the optical depth to re- 
scattering of cosmic-microwave-background (CMB) pho- 
tons suggests a reionization redshift z = 10.6 ± 1.4 [31 [5] 
if reionization occurred everywhere suddenly. However, 
the current paradigm is that reionization was highly in- 
homogeneous [B]: regions of ionized gas, seeded by the 
first light sources, grew and eventually coalesced, filling 
the entire volume of the IGM. If so, then reionization may 
have occurred over some redshift range. Recent searches 
for the redshifted 21-cm signal from neutral hydrogen 
[71 [8] now suggest that reionization occurred over a red- 
shift range Az > 0.06 [5]. The South Pole Telescope 
collaboration has now bounded that redshift range from 
above to be Az < 7.9 [TO], by searching for a fiuctu- 
ating kinetic-Sunyaev-Zeldovich (kSZ) signal fTT], under 
assumptions that the mechanism of reionization accords 
with prevailing theoretical models. 

Here we study the effects of patchy screening on the 
CMB dg. Thomson scattering of CMB photons during 
the EoR damps small-scale CMB fluctuations by a fac- 
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where T(n) is the optical depth in direction 
n on the sky. Patchy screening gives rise to a direction- 
dependent optical depth r(n). This then produces a B- 
mode polarization that is correlated in a characteristic 
way with the temperature and with the E mode polariza- 
tion [121 US] , and it also modulates the power in the tem- 
perature map. Here we interpret prior null searches for a 
modulation of CMB power in terms of an upper limit 
to optical-depth fluctuations, and we apply a minimum- 
variance estimator [T? for r fluctuations to the WMAP- 



7 temperature and polarization maps [T7| to search for 
patchy screening by measuring the off-diagonal TB cor- 
relations. We derive an upper limit to all multipoles of 
the power spectrum CJJ up to L = 512. We then discuss 
implications of these constraints for a simple phenomeno- 
logical reionization model whose parameters might serve 
as figures of merit for future experiments. We revisit pre- 
dictions for future experiments and discuss constraints on 
the parameter space imposed by the recent results from 
the EDGES |2j experiment. 

Formalism — Patchy screening suppresses primary 
anisotropics (marked with tilde), so the observed tem- 
perature fluctuation and polarization are, respectively, 



AT(n) = e-^(")Ar(n), 



p(n) = Q(n) + ^C/(n) 



P(n), 



(1) 



where Q and U are the usual Stokes parameters. All tem- 
perature and polarization correlations in the CMB can, 
in principle, be used to reconstruct the map of T(n). The 
EB estimator will ultimately provide the best sensitivity 
to patchy screening [T^l dH , once low-noise polarization 
measurements are available with future CMB exeriments. 
With WMAP and Planck |19j, however, the best sensitiv- 
ity is achieved with the TT correlation, which we discuss 
below. Here, we derive a constraint to patchy screening 
from the TB correlation, as a proof of principle. The 
estimator for the optical-depth fluctuation is [131 HI] > 
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where Yim and 2^im are spherical harmonics and spin- 
weighted spherical harmonics, respectively, and the sum 
is only over I + 1' + L = odd. The unbarred Bim 
and Tim are the observed temperature and polarization 
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multipoles, recovered from the maps and corrected for 
the combined instrumental-beam and pixeHzation trans- 
fer function Wi; bars represent the inverse- variance- 
filtered (IVF) multipoles, Bi„, = Bi,n/CP^ and fj. 



of fluctuations of r is 



Um 



Tim/Cf^ , where the TT and BB power spectra are ana- 
lytic estimates of the total (signal plus noise) power spec- 



trum Cf^^ 



-C; Wi m a given frequency 

band, for XX e {TT,BB}. The normalization can 
be calculated either analytically or using Monte Carlo 
simulations. The estimator is equivalent to the real part 
of the cosmic-birefringence estimator in Refs. [151 116] . 
the only difference being the parity condition. Ref. [T5] 
demonstrated that the full-sky formalism with the full- 
sky IVF procedure described above is justified in spite of 
the sky cuts introduced by masking the Galaxy. 
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Figure 1: The measurement of the power spectrum of fluctu- 
ations of the optical depth r from TB correlation with cor- 
responding 1(7 and 3(7 confidence-level intervals for all multi- 
poles up to the resolution limit of WMAP-7 is shown in the 
top two panels. A binned version with associated statisti- 
cal uncertainty is shown in the bottom panel. The first two 
bins are -0.0085 ± 0.1264 at L = 26, and 0.0029 ± 0.0056 at 
L — 77; they are omitted for the sake of clearer presentation. 
The measurements are consistent with zero at all multipoles. 



The TB correlations sought by this estimator can in 
principle also be generated by re-scattering of CMB pho- 
tons and by the kSZ effect from re-scattering. However, 
Ref. [T3] showed that the estimator is relatively insen- 
sitive to the kSZ effect, and also that only the large- 
scale {I <40) temperature fluctuations are sensitive to 
the former mechanism. In order to avoid large-scale con- 
tamination from pixel-pixel noise correlations in WMAP, 
we discard Tim and Eim multipoles below I = 100 from 
our analysis anyway, so we effectively probe only patchy 
screening. 

The estimator for the corresponding power spectrum 
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where /sky represents the fraction of the sky admitted 
by the analysis mask, correcting for the fact that the 
full-sky analysis is applied to the maps where a portion 
of the pixel values (mostly around the Galactic plane) 
was set to zero. When evaluated for the fixed cosmology 
of Ref. 4J and for the noise levels appropriate for the 
experiment in consideration, this four-point correlation 
provides a biased estimate of CJJ , where the bias mostly 
arises from the inhomogeneous pixel noise and the sky 
cuts. However, if this trispectrum is estimated by cross- 
correlating the tlm signal estimated from one frequency 
band with the same signal estimated from another fre- 
quency band, the largest contribution to its bias vanishes, 
because the instrumental noise is uncorrelated in differ- 
ent frequency bands. The leftover bias can be evaluated 
and subtracted by running a suite of null-hypothesis (no 
patchy signal) Monte Carlo simulations. We also use the 
null simulations to recover the statistical uncertainty for 
each measurement following the procedure described in 

Ref. nni. 
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Figure 2: Measurement of CJJ from TB correlation in 
WMAP-7 data. Results shown in black (filled circles) are 
obtained by using the analysis mask that covers all the point 
sources brighter than ~ 1 Jy, while the results in red (empty 
circles) are obtained after unmasking all the point sources. In 
spite of the large difference in the source contamination, the 
two results differ by much less than the statistical uncertainty, 
and no overall bias is observed. 

Results from TB estimator — We only show results for 
the cleanest band-cross-correlation [WW] [VV] where the 
estimate of tlm recovered from the W band is cross- 
correlated with the estimate from the V band. Prior to 
the analysis, we mask out the Galaxy and the known 
point sources using the fiducial 7-year analysis masks 
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Figure 3: Measurement of CJ^ from TB correlation in 
WMAP-7 data. Black filled circles represent the measure- 
ments obtained from the foreground-reduced maps after ap- 
plying the fiducial analysis mask (the fiducial result of Fig.[l|. 
The rest of the data points correspond to the two test cases: 
the green x's are obtained from the maps prior to foreground 
subtraction, but using the fiducial mask, while the red empty 
circles are measurements obtained from foreground-reduced 
maps after applying an extended mask. No overall bias is ob- 
served in the two cases, and all three results are consistent, 
within the estimated statistical uncertainty. 



available at the LAMBDA website [T^ (where for the 
combined mask /sky — 68%). After subtracting the bias, 
we recover a de-biased estimate CJJ of the power spec- 
trum at each multipole up to L = 512; Fig. [T] shows the 
binned measurements with estimated uncertainties. At 
all multipoles, we recover consistency with zero within 
the 3fT confidence level. 

Tests of Systematics — Our simulations do not include 
polarized point sources nor foreground residuals. In order 
to test their impact on our estimates of the power spec- 
trum and associated statistical uncertainty, we perform 
the tests described in Ref. [T5] . The results of these tests 
are shown in Figs. [2] and [3] which demonstrate that the 
foregrounds and point sources do not significantly affect 
the results reported in Fig. [ij 

Interpretation — To understand the implications of 
these measurements, we consider a simple parametriza- 
tion of inhomogeneous reionization in which optical- 
depth fluctuations are described by white noise smoothed 
on angular scales 9c = t^/Lc, or a power spectrum 



(47r/L^)(Ar)^e 



2-L''/Ll 



(4) 



shown in Fig. [4] for several values of Lc- We constrain 
the parameters At and Lc using the minimum-variance 
estimate [TS* for the amplitude. 



(A;? = (a[(Ar)2])2^C 



rr, fiducial 



Cr/var(Cr), (5) 



Figure 4: Family of simple patchy-reionization models, given 
by Eq. Q, for Ar = 1, and different values of Lc- 
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Figure 5: Shown are the la confidence-level upper limits from 
WMAP-7 TT and TB correlations on the amplitude [Arf of 
the patchy-screening model given in Eq. Q as a function of 
the coherence-scale parameter Lc- Also shown are sensitivity 
forecasts for experiments with map noise of 27^Karcmin and 
1/iKarcmin, and beam width of 7' and 5', corresponding to 
Planck and EPIC-like mission, respectively. The values ex- 
pected for a reionization surface that is crinkled on scales Rc 
are indicated by the gray band. Also shown is the portion of 
this parameter space excluded by EDGES [9]. 



where 



(a[(Ar)2 



«'^"""')Vvar(Cr) 
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is roughly the inverse-variance with which (Ar)^ can be 
measured, var(C£'^) is the variance of the power spec- 
trum, estimated from a suite of simulations with no 
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patchy screening, and CJ'^ are the unbinned measure- 
ments from WMAP-7 maps. Since the results are consis- 
tent with no signal, the variance provides a constraint on 
(Ar)^, which we show as a function of the model param- 
eter Lc in Fig. [5] We also show in Fig. [5] the upper limit 
to (Ar)^ inferred from upper limits to the power of XT 
modulation [14J. Given that the mean optical depth is 
known to be r ~ 0.1, it is clear that our bounds At < 1, 
from TB is far from constraining, and that At < 0.1, 
from TT is at best marginally constraining. 

Eq. ^ describes what happens if every point in the 
Universe goes suddenly from neutral to ionized, but with 
a reionization surface that is crinkled on a comoving scale 
of Rc ^ 200Mpc(Lc/150)"^ This smoothing scale, 
or bubble size, corresponds at a reionization redshift 

Zr ~ 10 to a redshift interval Az RcZr^^^lli'^ Hq/c. 
Since the optical depth scales with the reionization red- 
shift as T oc z^/^, we find that a bubble size Rc induces 
an optical-depth fluctuation At ~ 0.01(i?c/200Mpc). 
There is thus a rough scaUng, (At) ^ 0.01 {Lc/150)~^ , 
between the optical-depth-fluctuation amplitude and the 
correlation multipole Lc for the crinkly-surface model, 
represented by a thick band (to indicate roughly the the- 
ory uncertainty) in Fig. [5] Fig. [5] also shows the expec- 
tations [13l US] for the sensitivities of Planck and EPIC 
pUj . Also shown is a constraint for this crinkly-surface 
model from the EDGES constraint, Az > 0.06 (at 95% 
confidence) [9], from the all-sky redshifted 21-cm spec- 
trum. 

A wider range of reionization scenarios can be de- 
scribed by a "Swiss cheese" model in which bubbles of 
size Rc are spread over a larger redshift range iHj, so 
that each line of sight crosses, on average, N bubbles. 
The rms optical-depth fluctuation in Eq. Q would, for 
fixed Rc, then be reduced by a factor N^/^, relative to 
the crinkly-surface model. Thus, both the gray shaded 
area, and the "EDGES excluded" regions in Fig.[5]would 



be reduced by N^^"^. Note that kSZ fluctuations should 
increase in sensitivity as N increases [TOl [11] to comple- 
ment the reduced sensitivity of patchy screening in this 
limit. 

Conclusions — The directional dependence of the opti- 
cal depth T(n) encodes information about the morphol- 
ogy of the ionized regions during the epoch of reioniza- 
tion. Here we have used WMAP-7 temperature and po- 
larization data to derive a bound on the individual multi- 
poles of the optical-depth power spectrum up to i = 512, 
or bubble sizes larger than ~ 60 Mpc comoving. We then 
interpreted these null results in terms of a bound on an 
rms optical-depth fluctuation At in a model of white- 
noise fluctuations with coherence angle 6c- While the 
bound derived proves to be too weak to constrain realis- 
tic models, and probes bubble sizes larger than those 
(Rc ^ 10 Mpc) favored in current reionization mod- 
els, our result provides a proof of principle that such 
analyses can be carried out with future data. We then 
note that data from the forthcoming Planck satellite and 
from a subsequent post-Planck project should approach 
the realistic parameter space. Before such optical-depth- 
fluctuation searches are carried out in the future, though, 
several issues will need to be understood. For example, 
the estimator in Eq. ([2| has the same parity as that for 
the lensing potential [T5l [51] , and further modeling of the 
At and lensing signals, and/or de-lensing of the CMB, 
will be necessary for a CMB detection in the optical- 
depth fluctuation with Planck |19j or future-generation 
experiments. 
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